Primordial Nucleosynthesis 

By GARY STEIGMAN 

Department of Physics, The Ohio State University, Columbus, OH 43210, USA 

The primordial abundances of deuterium, helium-3, helium-4, and lithium-7 probe the baryon 
density of the Universe only a few minutes after the Big Bang. Of these relics from the early 
Universe, deuterium is the baryometer of choice. After reviewing the current observational sta- 
tus of the relic abundances (a moving target!), the baryon density determined by big bang 
nucleosynthesis (BBN) is derived. The temperature fluctuation spectrum of the cosmic back- 
ground radiation (GBR), established several hundred thousand years later, probes the baryon 
density at a completely different epoch in the evolution of the Universe. The excellent agreement 
between the BBN- and CBR-determined baryon densities provides impressive confirmation of 
the standard model of cosmology, permitting the study of extensions of the standard model. In 
combination with the BBN- and/or CBR-determined baryon density, the relic abundance of ''He 
provides an excellent chronometer, constraining those extensions of the standard model which 
lead to a nonstandard early-Universe expansion rate. 



1. Introduction 

As the hot, dense, early Universe rushed to expand and cool, it briefly passed through 
the epoch of big bang nucleosynthesis (BBN), leaving behind as relics the first com- 
plex nuclei: deuterium, helium-3, helium-4, and lithium-7. The abundances of these relic 
nuclides were determined by the competition between the relative densities of nucle- 
ons (baryons) and photons and, by the universal expansion rate. In particular, while 
deuterium is an excellent baryometer, "^He provides an accurate chronometer. Nearly 
400 thousand years later, when the cosmic background radiation (CBR) had cooled suffi- 
ciently to allow neutral atoms to form, releasing the CBR from the embrace of the ionized 
plasma of protons and electrons, the spectrum of temperature fluctuations imprinted on 
the CBR encoded the baryon and radiation densities, along with the universal expansion 
rate at that epoch. As a result, the relic abundances of the light nuclides and the CBR 
temperature fluctuation spectrum provide invaluable windows on the early evolution of 
the Universe along with sensitive probes of its particle content. 

The fruitful interplay between theory and data has been key to the enormous progress 
in cosmology in recent times. As new, more precise data became available, models have 
had to be reflned or rejected. It is anticipated this this process will - indeed, should - 
continue. Therefore, this review of the baryon content of the Universe as revealed by 
BBN and the CBR is but a signpost on the road to a more complete understanding of 
the history and evolution of the Universe. By highlighting the current successes of the 
present "standard" model along with some of the challenges to it, I hope to identify those 
areas of theoretical and observational work which will contribute to continuing progress 
in our endeavor to understand the Universe, its past, present, and future. 



2. A BBN Primer 

Discussion of BBN can begin when the Universe is already a few tenths of a second old 
and the temperature is a few MeV. At such early epochs the Universe is too hot and dense 
to permit the presence of complex nuclei in any significant abundances and the baryons 
(nucleons) are either neutrons or protons whose relative abundances are determined by 
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the weak interactions 



p + e < >n + i^e, 71 + < > p + ly^, n< >p + e + i^e- (2-1) 

The higher neutron mass favors protons relative to neutrons, ensuring proton dominance. 
When the weak interaction rates (Eq. 12. 1|) are fast compared to the universal expansion 
rate (and in the absence of a significant chemical potential for the electron neutrinos), 
n/p « exp(— Ato/T), where Am is the neutron-proton mass difference and T is the 
temperature {T^ — — T^, = Tn prior to annihilation). If there were an asym- 
metry between the number densities of Ve and ("neutrino degeneracy"), described 
by a chemical potential /ig (or, equivalently, by the dimensionless degeneracy parameter 

= /ie/r) then, early on, n/p « exp(— Am/T — ^e). For a significant positive chemical 
potential (^e ^ 0.01; more Ve than Pg) there are fewer neutrons than for the "standard" 
case (SBBN) which, as described below, leads to the formation of less ''He. 

The first step in building complex nuclei is the formation of deuterons via n+p < — >D+7. 
Sufhciently early on, when the Universe is very hot {T ^ 80 keV), the newly- formed 
deuterons find themselves bathed in a background of gamma rays (the photons whose 
relics have cooled today to form the CBR at a temperature of 2.7 K) and are quickly 
photo-dissociated, removing the platform necessary for building heavier nuclides. Only 
below ~ 80 keV has the Universe cooled sufficiently to permit BEN to begin, leading to 
the synthesis of the lightest nuclides D, ^He, "^He, and ^Li. Once BEN begins, D, ^H, 
and •^He are rapidly burned (for the baryon densities of interest) to "^He, the light nuclide 
with the largest binding energy. The absence of a stable mass-5 nuclide, in combination 
with Coulomb barriers, suppresses the BEN production of heavier nuclides; only ^Li is 
synthesized in an astrophysically interesting abundance. All the while the Universe is 
expanding and cooling. When the temperature has dropped below ~ 30 keV, at a time 
comparable to the neutron lifetime, the thermal energies of the colliding nuclides is insuf- 
ficient to overcome the Coulomb barriers, the remaining free neutrons decay, and BEN 
ends. 

From this brief overview of BEN it is clear that the relic abundances of the nuclides pro- 
duced during BEN depend on the competition between the nuclear and weak interaction 
rates (which depend on the baryon density) and the universal expansion rate (quanti- 
fied by the Hubble parameter H), so that the relic abundances provide early-Universe 
baryometers and chronometers. 



2.1. Early-Universe Expansion Rate 

The Friedman equation relates the expansion rate (measured by the Hubble parameter 
H) to the energy density (p): — where, during the early, "radiation-dominated" 

(RD) evolution the energy density is dominated by the relativistic particles present (p = 
Pr). For SEEN, prior to annihilation, these are: photons, pairs and, three flavors 
of left-handed (i.e., one helicity state) neutrinos (and their right-handed, antineutrinos). 

43 

Pr = Pj + Pe + 3p^ ^ —p^, (2.2) 

where p^ is the energy density in CBR photons. At this early epoch, when T ^ few MeV, 
the neutrinos are beginning to decouple from the 7 ~ plasma and the neutron to proton 
ratio, crucial for the production of primordial ^He, is decreasing. The time-temperature 
relation follows from the Friedman equation and the temperature dependence of p-y 



Pre - e=^ annihilation : t ^ 0.738 MeV^ s. 



(2.3) 



G. Steigman: Primordial Nucleosynthesis 



3 



To a very good (but not exact) approximation the neutrinos (ve, i^p., Vt) are decoupled 
when the pairs annihilate as the Universe cools below rrieC^. In this approximation 
the neutrinos don't share in the energy transferred from the annihilating pairs to the 
CBR photons so that in the post-e^ annihilation universe the photons are hotter than 
the neutrinos by a factor Tj/T^, = (ll/4)^/'^, and the relativistic energy density is 

PR — P-t + ipv — 1.68p^. (2.4) 

The post-e^ annihilation time-temperature relation is 

Post - e± annihilation : < = 1.32 MeV^ s. (2.5) 

2.1.1. Additional Relativistic Energy Density 

One of the most straightforward variations of the standard model of cosmology is to 
allow for an early (RD) nonstandard expansion rate H' = SH, where S = H' /H = t/t' 
is the expansion rate factor. One possibility for S" 7^ 1 is from the modification of the 
RD energy density (see Eas. 12.21 fc 12.41) due to "extra" relativistic particles X: pn 
PR + Px ■ If the extra energy density is normalized to that which would be contributed 
by one additional flavor of (decoupled) neutrinos (Steigman, Schramm & Gunn 1977), 
px = AN^p^ (N^ = 3 + AN^), then 

Spre = (.H'/H)pre = (1 + 0.163A7V,)i/2 ; Spost = [H' / H)post = (1 + 0.135A7V,)i/2. 

. ^.^-^^ 

Notice that S and lS.Ni, are related nonlinearly. It must be emphasized that it is S 
and not AA'^^ that is the fundamental parameter in the sense that any term in the 
Friedman equation which scales as radiation, decreasing with the fourth power of the 
scale factor, will change the standard-model expansion rate {S ^ 1). For example, higher- 
dimensional effects such as in the Randall-Sundrum model (Randall & Sundrum 1999a) 
may lead to either a speed-up in the expansion rate (5* > 1; AN^ > 0) or, to a slow- 
down (5' < 1; AN^ < 0); see, also, [Randall k, Sundrum (1999b)| [Binetruy et al. (2000)| 



Cline et al. (2000j 



2.2. The Baryon Density 

In the expanding Universe, the number densities of all particles decrease with time, so 
that the magnitude of the baryon density (or that of any other particle) has no meaning 
without also specifying when it is measured. To quantify the universal abundance of 
baryons, it is best to compare ne to the CBR photon density n-y. The ratio, ry = nB/n^ 
is very small, so that it is convenient to define a quantity of order unity, 

riio EE 10^°{nB/n,y) = 274f7B/i^ = 274 wb, (2.7) 

where JIb is the ratio (at present) of the baryon density to the critical density and h is 
the present value of the Hubble parameter in units of 100 km s~^ Mpc"'^ {cub = Ob^^)- 



3. BBN Abundances 

The relic abundances of D, ■^He, and ^Li are rate limited, determined by the competition 
between the early Universe expansion rate and the nucleon density. Any of these three 
nuclides is, therefore, a potential baryometer; see Figure ^ 

In contrast to the synthesis of the other light nuclides, once BBN begins (T 80 keV) 
the reactions building ''He are so rapid that its relic abundance is not rate limited. The 
primordial abundance of "^He is limited by the availability of neutrons. To a very good 
approximation, its relic abundance is set by the neutron abundance at the beginning 
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Figure 1. The SBBN-predicted abundances of D, ^He, and '^Li by number with respect to 
hydrogen, and the *He mass fraction Yp, as a function of the nucleon (baryon) abundance 
parameter 7710. The bands reflect the theoretical uncertainties (±lcr) in the BBN predictions. 

of BBN. As a result, the primordial mass fraction of ''He, Yp, while being a relatively 
insensitive baryometer (see Figure^, is an excellent, early-Universe chronometer. 

The qualitative effects of a nonstandard expansion rate on the relic abundances of the 
light nuclides may be understood with reference to Figure ^ For the baryon abundance 
range of interest the relic abundances of D and ^He are decreasing functions of 77; in this 
range, D and '^He are being destroyed to build ^He. A faster than standard expansion 
(5* > 1) provides less time for this destruction so that more D and ^He will survive. 
The same behavior occurs for ^Li at low values of 77, where its abundance is a decreasing 
function of ?/. However, at higher values of rj, the BBN-predicted ^Li abundance increases 
with 77, so that less time available results in less production and a smaller ^Li relic 
abundance. Except for dramatic changes to the early-Universe expansion rate, these 
effects on the relic abundances of D, "^He, and ^Li are subdominant to their variations with 
the baryon density. Not so for ^He, whose relic abundance is very weakly (logarithmically) 
dependent on the baryon density, but very strongly dependent on the early-Universe 
expansion rate. A faster expansion leaves more neutrons available to build ^He; to a 
good approximation AY ^ 0.16 {S — 1). It is clear then that if ^He is paired with any of 
the other light nuclides, together they can constrain the baryon density (rj or fie^^ = ^^b) 
and the early-Universe expansion rate {S or AN^). 

As noted above in |21 the neutron-proton ratio at BBN can also be modified from its 
standard value in the presence of a significant electron-neutrino asymmetry (^e ^ 0.01). 
As a result, Yp is also sensitive to any neutrino asymmetry. More Ve than De drives 
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the neutron-to-proton ratio down (see Ea. l2.1|l . leaving fewer neutrons available to build 
''He; to a good approximation AY « — 0.23^e (Kneller & Steigman 2003). In contrast, 
the relic abundances of D, ■^He, and ^Li are very insensitive to 7^ 0, so that when 
paired with ^He, they can simultaneously constrain the baryon density and the electron- 
neutrino asymmetry. Notice that if both S and are allowed to be free parameters, 
another observational constraint is needed to simultaneously constrain rj, S, and ^e- While 
neither '^He nor ^Li can provide the needed constraint, the CBR temperature anisotropy 
spectrum, which is sensitive to 77 and S, but not to ^g, can (see Barger et al. 2003b). 
This review will concentrate on combining constraints from the CBR and SBBN {S — 1, 
= 0) and also for 5" ^ 1 (^e = 0). For the influence of and constraints on electron 
neutrino asymmetry, see Barger et al. (2003b) and further references therein. 



4. Relic Abundances 

BBN constraints on the universal density of baryons and on the early-Universe expan- 
sion rate require reasonably accurate determinations of the relic abundances of the light 
nuclides. As already noted, D, "^He, and ''Li are all potential baryometers, while ^He is 
an excellent chronometer. The combination of the availablility of large telescopes and 
advances in detector technology has made it possible to obtain abundance estimates at 
various sites in the Galaxy and elsewhere in the Universe with unprecedented precision 
(statistically). However, the path to accurate primordial abundances is littered with sys- 
tematic uncertainties which have the potential to contaminate otherwise exquisite data. 
It is, therefore, fortunate that the relic nuclides follow very different post-BBN evolution- 
ary paths and are observed in diverse environments using a wide variety of astronomical 
techniques. Neutral deuterium is observed in absorption in the UV (or, in the optical 
when redshifted) against background, bright sources (O or B stars in the Galaxy, QSOs 
extragalactically) . Singly-ionized helium-3 is observed in emission in Galactic Hii re- 
gions via its spin- flip transition (the analog of the 21 cm line in neutral hydrogen). The 
helium-4 abundance is largely determined by observations of recombination lines of ion- 
ized (singly and doubly) ''He compared to those of ionized hydrogen in Galactic and, 
especially, extragalactic Hii regions. Observations of ^Li, at least those at low metallicity 
(nearly primordial) are restricted to absorption in the atmospheres of the oldest, most 
metal-poor stars in the halo of the Galaxy. The different evolutionary histories (described 
below) combined with the differrent observational strategies provide a measure of insur- 
ance that systematic errors in the determination of one of the light element abundances 
are unlikely to propagate into other abundance determinations. 

4.1. Deuterium - The Baryometer Of Choice 

The deuteron is the most weakly bound of the light nuclides. As a result, any deuterium 
cycled through stars is burned to "^He and beyond. Thus, its post-BBN evolution is 
straightforward: deuterium observed anywhere, anytime, provides a lower bound to the 
primordial D abundance. For "young" systems, in the sense of little stellar evolution (e.g. 
sites at high redshift and/or with very low metallicity), the observed D abundance should 
reach a plateau at the primordial value. Although there are observations of deuterium 
in the solar system and the interstellar medium (ISM) of the Galaxy which provide 
interesting lower bounds to its primordial abundance, the observations of relic D in a 
few, high redshift, low metallicity, QSO absorption line systems (QSOALS) are of most 
value in estimating its primordial abundance. 

While its simple post-BBN evolution is the greatest asset for relic D, the identical ab- 
sorption spectra of D i and H i (except for the velocity/wavelength shift resulting from the 
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Figure 2. Deuterium abundances, by number with respect to hydrogen D/H, versus metallicity 
(relative to solar on a log scale) from observations (as of early 2003) of QSOALS (filled circles). 
"X" is usually silicon or oxygen. Shown for comparison are the D abundances inferred for the 
local ISM (filled square) and the solar system (presolar nebula: "Sun" ; filled triangle) . 



heavier reduced mass of the deuterium atom) is a severe liability, limiting significantly 
the number of useful targets in the vast Lyman-alpha forest of the QSO absorption 
spectra (see Kirkman et al. (2003) for a discussion). It is essential in choosing a target 
QSOALS that its velocity structure be "simple" since a low column density H i absorber, 
shifted by ~ 81 km/s with respect to the main Hi absorber (an "interloper") would 
masquerade as D i absorption. If this is not recognized, a too high D /H ratio would be 
inferred. Since there are many more low-column density absorbers than those with high 
H I column densities, absorption systems with somewhat lower H i column density (e.g. 
Lyman-limit systems: LLS) are more susceptible to this contamination than the higher 
Hi column density absorbers (e.g. damped Lya absorbers: DLA). However, while the 
DLA have many advantages over the LLS, a precise determination of the Hi column 
density requires an accurate placement of the continuum, which could be compromised 
by interlopers. This might lead to an overestimate of the Hi column density and a con- 
comitant underestimate of D/H (J. Linsky, private communication). As a result of these 
complications, the path to primordial D using QSOALS has not been straightforward, 
and some abundance claims have had to be withdrawn or revised. Presently there are 
only five QSOALS with reasonably firm deuterium detections Kirkman et al. (2003) (and 
references therein); these are shown in Figure |2] along with the corresponding solar sys- 
tem and ISM D abundances. It is clear from Figure El that there is significant dispersion 
among the derived D abundances at low metallicity which, so far, mask the anticipated 
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Figure 3. Deuterium abundances versus the Hi column densities for the corresponding 

QSOALS shown in Figure H 



deuterium plateau. This suggests that systematic errors of the sort described here may 
have contaminated some of the determinations of the Di and/or Hi column densities. 

To explore the possibility that such systematic effects, which would be correlated with 
the Hi column density, may be responsible for at least some of the dispersion revealed 
in Figure 13 it is useful to plot the same QSOALS data versus the H i column density; 
this is shown in Figure |31 Indeed, there is the suggestion from this very limited data 
set that the low column density absorbers (LLS) have high D/H, while the high column 
density systems (DLA) have low abundances. However, on the basis of extant data it is 
impossible to decide which, if any, of these systems has been contaminated; there is no 
justification for excluding any of the present data. Indeed, perhaps the data is telling us 
that our ideas about post-BBN deuterium evolution need to be revised. 



To proceed further using the current data I follow the lead of O'Meara et al. (2001) 



and Kirkman et al. (2003) and adopt for the primordial D abundance the weighted mean 
of the D abundances for the five lines of sight (Kirkman et al. 2003); the dispersion in 
the data is used to set the error in j/d: 2/d = 2.6 ± 0.4. It should be noted that using the 
same data Kirkman et al. (2003) derive a slightly higher mean D abundance: ?/d = 2.74. 
The difference is traced to their first finding the mean of log(yD) and then using it to 
compute the mean D abundance (j/d = 10^'°^'^°'^). 

The BBN-predicted relic abundance of deuterium depends sensitively on the baryon 
density, t/d oc ?7~^'^, so that a ~ 10% determination of j/d can be used to estimate 
the baryon density to ~ 6%. For SBBN {S = I (N^ = 3), = 0), the adopted pri- 
mordial D abundance corresponds to ?7io(SBBN) = G-lOtoH i^Bh^ = 0.0223^0° olt), 
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in spectacular agreement with the 'Spergcl et al. (2003) estimate of 7710 — 6.14 ± 0.25 



{ilBh^ = 0.0224 ± 0.0009) based on WMAP and other GBR data (ACBAR and CBI) 
combined with large scale structure (2dFGRS) and Lynian-alpha forest constraints. In- 



deed, if the Spergel et al. (2003) estimate is used for the BBN baryon density, the BBN- 
predicted deuterium abundance is yu = 2.57 ± 0.27 (where a generous allowance of 
8% has been made for the uncertainty in the BBN prediction at fixed 77; for the 



Buries, Nollett & Turner (2001) nuclear cross sections and uncertainties the result is 

yu = 2mt'oi). 

4.2. Helium- 3 

Unlike D, the post-BBN evolution of ^He and ^Li are quite complex. '^He is destroyed in 
the hotter interiors of all but the least massive (coolest) stars, but it is preserved in the 
cooler, outer layers of most stars. In addition, hydrogen burning in low mass stars results 
in the production of significant amounts of new ^He (Iben 1967; Rood 1972; Dearborn, 
Steigman & Schramm 1986; Vassiliadis & Wood 1993; Dearborn, Steigman & Tosi 1996). 
To follow the post-BBN evolution of '^He, it is necessary to account for all these effects - 
quantitatively - in the material returned by stars to the interstellar medium (ISM). As 
indicated by the existing Galactic data (Geiss & Gloeckler 1998; Bania, Rood & Balser 
2002), a very delicate balance exists between net production and net destruction of ^He in 
the course of the evolution of the Galaxy. As a consequence, aside from noting an excellent 
qualitative agreement between the SBBN predicted and observed '^He abundances, ^He 
has - at present - little role to play as a quantitatively useful baryometer. In this spirit, it 
is noted that an uncertain estimate of the primordial abundance of '^He may be inferred 
from the observation of an outer-Galaxy (less evolved) Hii region (Bania et al. 2002): 
j/3 = 105(3He/H) = 1.1 ±0.2. 

4.3. Lithium-7 

A similar scenario may be sketched for ^Li. As a weakly bound nuclide, it is easily 
destroyed when cycled through stars except if it can be kept in the cooler, outer layers. 
The high lithium abundances observed in the few "super-lithium-rich red giants" provide 
direct evidence that at least some stars can synthesize post-BBN lithium and bring it 
to the surface. But, an unsolved issue is how much of this newly-synthesized lithium is 
actually returned to the ISM rather than mixed back into the interior and destroyed. 

With these caveats in mind, in Figure 01 lithium abundances are shown as a function 
of metallicity from a compilation by V. V. Smith (private communication). Since the 
quest for nearly primordial lithium is restricted to the oldest, most metal-poor stars in 
the Galaxy, stars that have had the most time to redistribute - and destroy or dilute - 
their surface lithium abundances, it is unclear whether the "plateau" at low metallicities 
is representative of the primordial abundance of lithium. Although it seems clear that 
the lithium abundance in the Galaxy has increased since BBN, a quantitatively reliable 
estimate of its primordial abundance eludes us at present. Given this state of affairs, the 
most fruitful approach is to learn about stellar structure and evolution by comparing the 
BBN-predicted lithium abundance to those abundances inferred from observations of the 
oldest stars, rather than to attempt to use the stellar observations to constrain the BBN- 
inferred baryon density. Goncentrating on the low-metallicity, nearly primordial data, 
it seems that [Li] = 12-|-log(Li/H) « 2.2 ± 0.1. This estimate will be compared to the 
BBN-predicted lithium abundance using D as a baryometer and, to the BBN-predicted 
lithium abundance using the GBR-inferred baryon density. Any tension between these 
BBN-predicted abundances and that inferred from the Galactic data may provide hints 
of nonstandard stellar astrophysics. 
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• McDonald Spectra analyzed here 
□ Lambert et al. (1991) 
A Cunhaetal. (1995) 
O NIssen etal. (1999) 
o Ryan et al. (1998) 
o Ryan et al. (1999) 
SS= Solar System 
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Figure 4. Lithium abundances, log e(Li) = [Li] = 12+log(Li/H) versus metallicity (on a log 
scale relative to solar) from a compilation of stellar observations by V. V. Smith. 



4.4. Helium-4 - The BBN Chronometer 

The good news about ^He is that, as the second most abundant nuclide, it may be 
observed throughout the Universe. The bad news is that its abundance has evolved 
since the end of BBN. In order to infer its primordial value it is therefore necessary to 
track the ^He abundance determinations (mass fraction Yp) as a function of metallicity 
or, to limit observations to very low metallicity objects. Although, as for D, there are 
observations of ''He in the ISM and the solar system, the key data for determining its 
primordial abundance comes from observations of metal-poor, extragalactic Hii regions. 
A compilation of current data (courtesy of K. A. Olive) is shown in Figure [S] where 
the ^He mass fraction is plotted as a function of the oxygen abundance; note that the 
solar oxygen abundance, O/H « 5 x 10~* (AUende-Prieto, Lambert & Asplund 2001) is 
off-scale in this figure. These are truly low metallicity H ii regions. 

It is clear from Figure|31that the data exist to permit the derivation of a reasonably ac- 
curate estimate (statistically) of the primordial ''He mass fraction Yp, with or without any 
extrapolation to zero-metallicity. What is not easily seen in Figure[Slgiven the Yp scale, is 



that Yp derived from the data assembled from the literature by Olive & Steigman (1995) 



and Olive, Skillman & Steigman (1997) (Yp = 0.234 ± 0.003) is marginally inconsistent 



(at ~ 2a) with the value derived by Izotov, Thuan fc Lipovetsky (1997) and Izotov fc Thuan (1998) 
from their nearly independent data set (Yp = 0.244 ± 0.002). In addition, there are a 
variety of systematic corrections which might modify both data sets (Steigman, Viegas & 
Gruenwald 1997; Viegas, Gruenwald & Steigman 2000; Olive & Skillman 2001; Sauer & 
Jedamzik 2002; Gruenwald, Steigman & Viegas 2002; Peimbert, Peimbert & Luridiana 
2002) 

Unless/until the differences in Yp derived by different authors from somewhat dif- 
ferent data sets is resolved and the known systematic errors are corrected for (the un- 
known ones will always hang over us like the sword of Damocles), the following com- 
promise, adopted by Olive, Steigman fc Walker (2000)[ may not be unreasonable. From 
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Figure 5. The ''He mass fraction Y derived from observations of extragalactic H ii regions of 
low metallicity versus the corresponding Hii region oxygen abundances (from K. A. Olive). 



Olive & Steigman (1995) and Olive, Skillman & Steigman (1997)1 t he 2a range for Yp is 
0.228 - 0.240, while from the |Izotov, Thuan fc Lipovetsky (1997)| and |!^ov fc Thuan (1998) 
data the 2a range is Yp = 0.240 - 0.248. Thus, although the current estimates are likely 
dominated by systematic errors, they span a ~ 2cr range from Yp = 0.228 to Yp = 
0.248. Therefore, as proposed by Olive, Steigma n fc Walker (2000)[ we adopt here a cen- 
tral value for Yp = 0.238 and a ~ Ict uncertainty of 0.005: Yp = 0.238 ± 0.005. Given 
the approximation (see ^ AY « 0.16 (5 — 1), for ctyp ~ 0.005 the uncertainty in S is 
~ 0.03 (corresponding to an uncertainty in AA''^ of « 0.4). 



5. The Baryon Density From SBBN 

For SBBN, where S = I {N^ — 3) and = 0, the primordial abundances of D, 
^He, ^He, and ^Li are predicted as a function of only one free parameter, the baryon 
density parameter {q or ^l-^h? = wb). As described above (see ii4.1() . D is the baryometer 
of choice. From SBBN and the adopted relic abundance of deuterium, j/d = 2.6 ± 0.4, 
= e.ltas {^Bh"^ = 0.022 ± 0.002). 

Having determined the baryon density to ~ 10% using D as the SBBN baryometer, 
it is incumbent upon us to compare the SBBN-predicted abundances of the other light 
nuclides with their relic abundances inferred from the observational data. For this baryon 
density, the predicted primordial abundance of '^He is — 1.04 ± 0.10, in excellent 
agreement with the primordial value of ys = 1.1 ± 0.2 inferred from observations of an 
outer-Galaxy Hii region (Bania et al. 2002). Within the context of SBBN, D and "^He 
are completely consistent. 

The first challenge to SBBN comes from ''He. For the SBBN-determined baryon density 
the predicted ''He primordial mass fraction is Yp = 0.248 ± 0.001, to be compared with 
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our adopted value from extragalactic Hii regions (Olive, Steigman & Walker 2000) of 
Y^sw ^ 0.238 ± 0.005. Agreement is only at the ~ 2a level. Given the unresolved 
systematic uncertainties in determining Yp from the Hii region data, it is not clear 
at present whether this is a challenge to SBBN or to our understanding of H ii region 
recombination spectra. As will be seen below, this tension between SBBN D and "^He can 
be relieved for nonstandard BBN if the assumption that S — 1 {N^ ~ 3) is relaxed. 

The conflict with the inferred primordial abundance of lithium is even more challenging 
to SBBN. For = 2.6 ± 0.4, [Li] = 2.65lo ;[2. This is to be compared to the estimate 
(see Figure of [Li] = 2.2 ± 0.1 based on a sample of metal-poor, halo stars. The 
conflict is even greater with the Ryan et al. (2000) estimate of [Li] = 2.G9j]Q |^3 derived 
from an especially selected data set. Unlike the tension between SBBN and the D and 
^He abundances, the conflict between D and ^Li cannot be resolved by a nonstandard 
expansion rate (nor, by an electron neutrino asymmetry). Most likely, the resolution 
of this conflict is astrophysical since the metal-poor halo stars from which the relic 
abundance of lithium is inferred have had the longest time to mix their surface material 
with that in their hotter interiors, diluting or destroying their prestellar quota of lithium 



(see, e.g. Pinsonneault et al. (2002) and references to related work therein). 

At present SBBN in combination with the limited data set of QSOALS deuterium 
abundances yields a ~ 10% determination of the baryon density parameter. Consistency 
between the inferred primordial abundances of D and '^He lends support to the internal 
consistency of SBBN, but the derived primordial abundances of ^He and ^Li pose some 
challenges. For ^He the disagreement is only at the ~ 2a level and the errors in the obser- 
vationally inferred value of Yp are dominated by poorly quantifled systematics. However, 
if the current discrepancy is real, it might be providing a hint at new physics beyond the 
standard model (e.g. nonstandard expansion rate and/or nonstandard neutrino physics). 
Before considering the effects on BBN of a nonstandard expansion rate (S* ^ 1; N^, 7^ 3), 
we will compare the SBBN estimate of the baryon density parameter with that from the 
GBR. 



6. The Baryon Density From The CBR 

Some 400 kyr after BBN has ended, when the Universe has expanded and cooled suf- 
ficiently so that the ionized plasma of protons, alphas, and electrons combines to form 
neutral hydrogen and helium, the CBR photons are set free to propagate throughout the 
Universe. Observations of the CBR today reveal the anisotropy spectrum of temperature 
fluctuations imprinted at that early epoch. The so-called acoustic peaks in the tempera- 
ture anisotropy spectrum arise from the competition between the gravitational potential 
and the pressure gradients. An increase in the baryon density increases the inertia of 
the baryon - photon fluid shifting the locations and the relative heights of the acoustic 
peaks. In Figure |^ are shown three sets of temperature anisotropy spectra for cosmo- 
logical models which differ only in the choice of the baryon density parameter ujb- Also 
shown in FigureElare the WMAP data from Bennett et al. (2003) It is clear from Figure 
El that the CBR provides a very good baryometer - independent of that from SBBN and 
primordial deuterium. Based on the WMAP data alone, Barger et al. (2003a) find that 
the best fit value for the density parameter is 7710 = 6.3 (wb = 0.023) and that the 2a 
range extends from 7710 — 5.6 to 7.3 (0.020 ^ wb =^ 0.026). This is in excellent (essentially 



perfect!) agreement (as it should be) with the CBR-only result of Spergel et al. (2003) 



More importantly, as may be seen clearly in Figured (courtesy of D. Marfatia), this 
independent constraint on the baryon density parameter, sampled some 400 kyr after 
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Figure 6. The CBR temperature fluctuation anisotropy spectra for three choices of the baryon 
density parameter ujb ~ 0.018, 0.023, 0.028, in order of increasing height of the first peak. Also 
shown are the WMAP data points. 
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Figure 7. The normalized likelihood distributions for the baryon density parameter 7710 derived 
from SBBN and the primordial abundance of deuterium (solid curve; see H4.H and from the CBR 
using WMAP data alone (dashed curve). The bottom horizontal axis is the baryon-to- photon 
ratio parameter 7710; the top axis is the baryon density parameter uib = flsh^ ■ 
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Figure 8. The BBN-predicted relation between the *He mass fraction Yp and the deuterium 
abundance yo for three, early-Universe expansion rates corresponding to A'^^ — 2, 3, 4. The filled 
circle with error bars is for the D and *He primordial abundances adopted here. 

BBN, is in excellent agreement with that from SBBN (see ^S)), providing strong support 
for the standard model of cosmology. 

The independent determination of the baryon density parameter by the CBR rein- 
forces the tension between SBBN and the relic abundances of '^He and ^Li inferred from 
the observational data (see fjS)). In the context of SBBN, the slightly higher best value of 
r] from the WMAP data (compared to that from D plus SBBN) increases the expected 
primordial abundances of ^He and ""Li (see Figure^, widening the gaps between the 
SBBN predictions and the data. Keeping in mind the observational and theoretical diffi- 
culties in deriving the primordial abundances from the data, it is nonetheless worthwhile 
to explore a class of nonstandard alternatives to the standard model of cosmology in 
which the early Universe expansion rate is modified (5 7^ 1, N^, ^ 3). 

7. Nonstandard BBN: S" / 1, iV^ / 3 

As outlined in ^ for fixed rj a.s S increases the BBN-predicted abundances of D, ■^He, 
and ''He increase (less time to destroy D and ^He, more neutrons available for ''He), 
while that of ^Li decreases (less time to produce ^Li). Since it is the "^He mass fraction 
that is most sensitive to changes in the early Universe expansion rate and, since the 
SBBN-predicted value of Yp is too large when compared to the data, S < 1 {Ni, < 3) is 
required. For a slower than standard expansion rate the predicted abundances of D and 
^He decrease compared to their SBBN values (at fixed rj) while that of ^Li increases. Since 
the BBN-predicted abundance of D increases with decreasing baryon density, a decrease 
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Figure 9. The la, 2a, and 3a contours in the AN^, - rjio plane from BBN and the rehc D and 
*He abundances. The best fit values of AA*"^ and ryio are marked by the cross. 



in S can be compensated for by a decrease in 77. For 7710 « 6 and S' — 1 <C 1, a good 
approximation (for fixed D) is ArjiQ « 6(5' — 1) (Kneller & Steigman 2003). In Figure 
|H1 are shown the ^He - D (Yp versus D /H) relations for three values of the expansion 
rate parameterized by N^. To first order, the combination of 77 and S that recovers the 
SBBN deuterium abundance will leave the '^He abundance prediction unchanged as well, 
preserving its good agreement with the observational data. However, the consequences 
for ^Li are not so favorable. The BBN abundance of ^Li increases with decreasing S 
but decreases with a smaller 77; the two effects nearly cancel leaving essentially the same 
discrepancy as for SBBN. For ''Li, a nonstandard expansion rate cannot relieve the tension 
between the BBN prediction and the observational data. 

Setting aside ^Li, it is of interest to consider the simultaneous constraints from BBN 
on the baryon density parameter and the expansion rate factor from the abundances of 
D and ^He; it has already been noted that for this nonstandard case, D and '^He will 
remain consistent. In Figure |H| are shown the la, 2a, and 3a contours in the AA^^, - 77 
plane derived from BBN and the D and "^He relic abundances. As expected from the 
discussion above, the best fit value of 77 (the cross in Figure IHJl has shifted downward to 
7710 = 5.7 {lob = 0.021). While the best fit is for AN^ = -0.7 {S = 0.94), it should be 
noted that the standard case of N^, = 3 is entirely compatible with the data at the ^ 2a 
level. 



8. Nonstandard CBR: S ^ I, ^ 3 

The CBR temperature fluctuation anisotropy spectrum is sensitive to the early-Universe 
radiation density (pr) as well as to the overall expansion rate. The early Universe is ra- 
diation dominated so that p « pR oc 1 -I- 0.135AN^ (see Eg. 12. 61 and recall that p oc H"^). 
The late Universe is matter dominated (MD) {ujm = ^uh^) and the crossover from RD 
to MD, important for the growth of fluctuations and for the age/size of the Universe at 
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recombination, occurs for a redshift 

zeq = 2.4 X W^LOuil + 0.135AN^)-^ (8.1) 

If the matter content is kept fixed while the radiation content is increased, corresponding 
to a faster than standard expansion rate, matter-radiation equahty is delayed, modifying 
the growth of fluctuations prior to recombination and, also, the Universe is younger 
at recombination and has a smaller sound horizon, shifting the angular location of the 
acoustic peaks. The degeneracy between the radiation density (AA^^ or S) and u>m is 
broken by the requirement that Om + = 1 and the HST Key Project determination 



of the Hubble parameter (see Barger et al. (2003a) for details and further references). 



In Figure El are shown the la and 2a contours in the AN^ - rj plane from the CBR 
(WMAP) data; note the very different AN^, scales and ranges in Figures El and EH As 
is the case for BEN (see fJZjl, the CBR favors a slightly slower than standard expansion. 
However, while the "best" fit value for the expansion rate factor is at S' < 1 {AN^, < 0), 
the CBR likelihood distribution of AA^;^ values is very shallow and the WMAP data are 
fully consistent with 5=1 (AA^^ = 0). 

Comparing Figures li^ and [TOl it is clear that for this variant of the standard cosmology 
there is excellent overlap between the rj - ANi, confidence contours from BBN and those 
from the CBR (see Barger et al. 2003a). This variant of SBBN {S ^ 1) is consistent with 
the CBR. In Figure ^2 (from Barger et al. 2003a) are shown the confidence contours 
in the 77 - AN^, plane for a joint BBN - CBR fit. Again, while the best fit value for 
AA'^i^ is negative (driven largely by the adopted value for Yp), AA^^ = (S" = 1) is quite 
acceptable. 
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Figure 11. As for Figure ITI71 but for the joint BBN - CBR fit. The best fit point (r;io = 6.0, 
ANi, — —0.75) is indicated by the cross. 



9. Summary And Conclusions 

As cosmology deals with an abundance of precision data, redundancy will be the key 
to distinguishing systematic errors from evidence for new physics. BBN and the CBR 
provide complementary probes of the Universe at two epochs widely separated from each 
other and from the present. For the standard model assumptions (N^, — 3, S = 1) the 
SBBN-inferred baryon density is in excellent agreement with that derived from the CBR 
(with or without the extra constraints imposed by large scale structure considerations 
and/or the Lyman alpha forest). For this baryon density (7710 « 6.1, wb ~ 0.022), the 
SBBN-predicted abundances of D and '^He are in excellent agreement with the observa- 
tional data. For ''He the predicted relic mass fraction is ~ 2a higher than the primordial 
abundance inferred from current data, hinting at either new physics or the presence of 
unidentified systematic errors. For ^Li too, the SBBN-predicted abundance is high com- 
pared to that derived from very metal-poor stars in the Galaxy. While the tension with 
*He can be relieved by invoking new physics in the form of a nonstandard (slower than 
expected) early-Universe expansion rate, this choice will not reconcile the BBN-predicted 
and observed abundances of ^Li. When both the baryon density and the expansion rate 
factor are allowed to be free parameters, BBN (D, ^He, and ''He) and the CBR (WMAP) 
agree at 95% confidence for 5.5 ^ 7710 ^ 6.8 and 1.65 ^ Nj^ < 3.03. 

The engine powering the transformation of the study of cosmology from its youth to 
its current maturity has been the wealth of observational data accumulated in recent 
years. In this data-rich, precision era BBN, one of the pillars of modern cosmology, con- 
tinues to play a key role. The spectacular agreement between the estimates of the baryon 
density derived from processes at widely separated epochs has confirmed the general 
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assumptions of the standard models of cosmology and of particle physics. The tension 
with "^He (and with ^Li) provides a challenge, along with opportunities, to cosmology, 
to astrophysics, and to particle physics. Whether the resolution of these challenges is 
observational, theoretical or, a combination of both, the future is bright. 
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presented here are from recent collaborations or discusions with V. Barger, J. P. Kneller, 
H.-S. Lee, J. Linsky, D. Marfatia, K. A. Olive, R. J. Scherrer, V. V. Smith, and T. 
P. Walker. My research is supported at OSU by the DOE through grant DE-FG02- 
91ER40690. This manuscript was prepared while I was visiting the Instituto Astronomico 
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